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A  small  thermoelectric  generator  to  power  autonomous  sensors  in  remote  environmental  sites  is  studied, 
designed,  realized,  characterized,  and  tested.  The  thermoelectric  phenomena,  applied  to  our  device,  are 
theoretically  introduced  and  experimentally  verified  by  directly  measuring  the  physical  quantities  when 
the  thermoelectric  generator  operates  in  working  conditions.  The  device  is  then  tested  under  different 
external  conditions,  showing  that  it  is  able  to  supply,  for  sufficient  long  time,  an  output  voltage  higher 
than  200  mV  and  an  output  power  on  the  order  of  10  mW  when  a  temperature  difference  higher  than 
10  K  and  a  load  resistance  close  to  the  internal  resistance  are  considered.  Furthermore  we  developed  a 
devoted  power  conditioning  circuit  in  order  to  usefully  manage  the  output  voltage.  Finally,  we  tested 
the  device  in  real  operative  conditions. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

During  the  last  20  years  the  request  of  renewable  and  environ¬ 
mentally  “green”  sources  of  energy  in  order  to  deal  with  the  prob¬ 
lem  of  pollution  and  global  warming  has  been  noticeably 
increased.  In  particular,  parts  of  these  studies  are  focused  on  the 
pursuit  of  energy  harvesting  from  unused  sources  like  mechanical 
vibrations  and  wasted  heat.  Within  these  researches,  the  applica¬ 
tions  which  involve  the  thermoelectric  phenomena  have  been 
intensively  developed  [1-5].  In  more  details,  the  thermoelectric 
modules  are  solid-state  devices  commonly  used  in  cooling/heating 
applications,  but  in  the  last  years  they  have  been  largely  studied 
and  utilized  as  energy  converters  since  they  are  able  to  produce 
electricity  from  heat  absorption.  In  particular,  small  thermoelectric 
generator  devices  are  principally  used  to  power  localized  autono¬ 
mous  sensors,  but  their  fields  of  application  are  at  this  time  re¬ 
stricted  to  aero-spatial  and  automotive  industries. 

Within  this  framework,  in  this  paper  we  extend  the  area  of 
interest  for  the  employment  of  thermoelectric  generators  to  the 
environmental  monitoring.  The  use  of  large  thermoelectric  gener¬ 
ators  in  remote  sites  started  since  the  sixties  of  the  previous  cen¬ 
tury  in  order  to  power  data  acquisition  equipments  and 
communication  devices,  such  as  weather  stations,  navigation  aids, 
and  subsea  operations  [6];  in  all  these  cases  the  heat  source  was 
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powered  by  either  radioisotopes  or  fossil  fuels.  In  our  work  we 
studied,  designed,  realized,  and  characterized  a  small  thermoelec¬ 
tric  generator  in  order  to  supply  energy  to  autonomous  sensors, 
placed  in  remote  locations  with  hard  environmental  conditions 
(very  low  temperature,  difficult  access,  etc.),  where  the  usual 
renewable  energy  sources,  like  solar  and  wind  energy,  are  not  reg¬ 
ularly  available. 

The  novelty  of  our  work,  with  respect  to  the  previously  men¬ 
tioned  devices  and  the  known  applications  discussed  in  literature, 
has  to  be  found  both  in  the  originality  of  the  proposed  device,  by 
considering  in  particular  its  whole  structure,  and  the  field  of  appli¬ 
cation.  Furthermore,  this  thermoelectric  generator  has  been  com¬ 
pletely  realized,  optimized,  and  characterized  after  having 
considered  both  the  physical  models  for  thermoelectric  phenom¬ 
ena  and  its  specific  practical  applications.  As  first  consideration, 
we  highlight  the  peculiarity  of  the  heat  source  employed  since  in 
our  case  the  heat  is  supplied  by  a  flameless  catalytic  burner.  Such 
type  of  heat  source  has  been  already  studied  [7,10],  but  only  in  lab¬ 
oratory  tests,  without  any  sort  of  applications  in  real  operative 
conditions.  This  kind  of  heat  source  is  more  competitive  and  reli¬ 
able  with  respect  to  flame  burner  when  it  is  used  for  applications 
where  low  power  and  small  dimensions  are  required.  In  addition, 
our  study  permits  to  model  our  thermoelectric  generator  when  dif¬ 
ferent  natural  heat  sources  are  considered  since  the  employed  cat¬ 
alytic  burner  basically  provides  a  low-temperature  heat  source, 
exactly  what  we  expect  from  renewable  heat  sources  as  geother¬ 
mal,  ocean  thermal  energy,  and  waste  heat  to  mention  few  ones 
[6].  Secondly,  the  comparison  between  the  physical  model  and 
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the  experimental  studies  [15,14,11,9]  of  thermoelectric  generators 
as  well  as  their  performances  [8]  under  different  operative  condi¬ 
tions  have  been  extensively  developed  but  only  in  the  ideal  labora¬ 
tory  conditions.  In  our  work,  all  the  experimental  tests  have  been 
worked  out  by  using  the  device  in  its  real  operative  layout.  More¬ 
over,  the  situations  actually  found  in  realistic  operative  conditions 
are  mainly  studied  theoretically  [12],  whereas  the  realistic  experi¬ 
mental  tests  are  principally  developed  by  considering  residential 
applications  [13]  or  devices  assembled  with  large  parts  [6].  Here, 
we  have  extensively  tested  our  small  prototype  under  realistic 
environmental  conditions,  the  same  ones  which  it  can  find  when 
it  effectively  operates  in  the  required  applications.  Finally,  we  re¬ 
mark  that  the  power  conditioning  electric  module  is  generally  ne¬ 
glected  in  thermoelectric  generator  demonstrations,  whereas  in 
this  work  it  plays  a  fundamental  role. 

This  paper  is  organized  as  it  follows.  After  having  introduced 
the  thermoelectric  phenomena  we  recall  their  basic  physical  mod¬ 
el.  Subsequently,  we  describe  the  design  of  the  whole  thermoelec¬ 
tric  generator.  In  particular,  we  discuss  in  details  the  power 
conditioning  external  circuit,  developed  in  order  to  correctly  use 
the  output  voltage  as  energy  supply  for  autonomous  sensors.  Then, 
we  illustrate  the  experimental  measurements  carried  out  both  to 
test  the  validity  of  the  physical  model  previously  discussed  and 
to  check  the  capability  of  our  device  to  be  reasonably  used  in  the 
expected  environment.  Finally,  we  show  the  results  obtained  dur¬ 
ing  the  tests  carried  out  in  real  operative  conditions.  By  summariz¬ 
ing,  we  have  designed  and  realized  a  complete  thermoelectric 
generator  powered  by  a  flameless  catalytic  burner;  the  device 
has  been  tested  both  in  laboratory,  by  comparing  the  experimental 
results  with  the  theoretical  predictions,  and  in  real  environmental 
situations,  by  proving  its  capability  to  be  used  as  power  generator 
for  environmental  monitoring  sensors. 


2.  Physical  models 

Thermoelectric  (TE)  phenomena  are  basically  observed  when 
two  different  materials  are  joined  together  and  they  are  either  sub¬ 
jected  to  a  temperature  gradient  or  crossed  by  an  electrical  current 
[6].  In  particular,  the  Seebeck  effect  describes  the  generation  of  a 
voltage  difference  AV  when  the  two  sides  of  the  junction  are  held 
at  different  temperature  AT,  showing  a  proportional  dependence 
given  by  AV  =  aAT,  where  a  is  the  Seebeck  coefficient  or  thermo¬ 
power.  The  inverse  effect,  the  Peltier  effect,  arises  when  an  electri¬ 
cal  current  passes  through  the  junction  of  two  different  materials, 
with  the  consequent  heat  absorption  or  rejection,  depending  on 
the  current  direction.  The  heat  flow  is  proportional  to  the  current 
7  as  dQJdt  =  77/,  being  77  the  Peltier  coefficient  which  Thomson  (Lord 
Kelvin)  proved  to  be  proportional  to  the  absolute  temperature 
77  =  ocT.  A  further  TE  effect,  the  Thomson  effect,  describes  the  heat 
production  (absorption  or  emission)  of  an  electrical  current  in  a 
single  material  under  the  presence  of  a  temperature  gradient.  By 
the  way,  the  Seebeck,  Peltier,  and  Thomson  effects  are  reversible, 
whereas  the  Joule  effect  is  obviously  not. 

Modern  TE  refrigeration  devices  and  TE  power  generation  de¬ 
vices  are  based  on  Peltier  and  Seebeck  effects,  respectively;  their 
core  is  made  by  thermoelectric  couples  composed  of  n-type  and 
p-type  semiconductors  (basically,  Bi2Te3  and  Si^Ge*)  pellets  con¬ 
nected  by  metallic  electrical  contact  pads.  In  such  devices,  a  large 
number  of  these  thermoelectric  couples  are  connected  electrically 
in  series  and  thermally  in  parallel.  Thermocouples  are  then  placed 
between  two  ceramic  plates  which  are  electrically  insulating  and 
thermally  conducting.  An  example  of  TE  module  (TEM)  is  shown 
in  Fig.  1.  Thermal  energy  is  supplied  by  the  heat  source  at  temper¬ 
ature  Th  which  drops  at  temperature  Thj  at  the  hot  thermocouples 
junctions  due  to  the  finite  thermal  conductance  I<  of  the  ceramic 


Th 

r~ 


Tc 


Fig.  1.  Sketch  of  TEM.  Grey  surfaces:  ceramic  plates;  +/-:  electrode  terminals; 
yellow  lines:  electrical  connections;  black  and  red  bars:  n-  and  p-type  semicon¬ 
ductors;  Th :  temperature  of  heat  source;  Tc:  temperature  of  heat  sink;  W- 
temperature  of  hot  junction;  Tcj:  temperature  of  cold  junction.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  this  article.) 

plates.  For  the  same  reason,  the  cold  junction  remains  at  the  tem¬ 
perature  TCJ,  whereas  at  the  heat  sink  we  recover  the  lower  tem¬ 
perature  Tc.  We  therefore  define  the  two  temperature  differences 
[14]  as  AT  =  Tjj  —  Tc  and  ATeff  =  T^j  —  Tcj . 

Geometrically,  TEMs  are  characterized  by  the  cross  sectional 
area  A  of  ceramic  plates,  the  thickness  77  of  the  whole  TEM,  the 
cross  sectional  area  As  and  the  height  h  of  each  pellet,  the  height 
h'  =  (77  -  h)l 2  of  each  ceramic  plate,  and  the  number  N  of  thermo¬ 
couples.  The  thermal  conductance  I<  of  ceramic  plates  is  given  by 
I<  =  kcA/h\  where  kc  is  the  thermal  conductivity  of  the  ceramic 
whose  value  is  in  the  range  (18-36)  W/Km  for  A1203,  typically  used 
in  TEMs  applications.  Besides,  the  other  properties  of  the  TEM  are 
the  internal  thermal  conductance  I<in  =  2 NksAslh  and  the  internal 
electrical  resistance  Rin  =  2NpshlAs ,  where  ks  (typically  1. 5-2.0  W / 
Km)  and  ps  are  the  thermal  conductivity  and  the  electrical  resistiv¬ 
ity  of  a  single  pellet,  respectively  (for  the  corresponding  values  in 
literature  see,  e.g.,  [6,16,15]).  Without  loss  of  generality,  we  have 
considered  the  same  thermal  and  electrical  properties  for  both  n- 
and  p-type  semiconductors;  different  properties  can  be  eventually 
taken  into  account. 

In  the  sequel,  we  will  consider  the  simplest  and  most  com¬ 
monly  used  model  to  test  an  ideal  TEG  (see  Ref.  [6]  for  a  detailed 
discussion).  This  classical  method  is  based  on  the  assumption  that 
all  the  physical  coefficients  are  constant  and  temperature  indepen¬ 
dent,  consequently  the  Thomson  effect  is  neglected.  Although  its 
simplicity,  this  method  has  been  successfully  used  to  model  a  large 
variety  of  applications  (e.g.,  see  Refs.  [15,14]).  Nevertheless,  we  no¬ 
tice  that  more  complete  and  detailed  models  have  been  exten¬ 
sively  studied  and  experimentally  verified  [16,17].  On  the  other 
hand,  this  analysis  is  beyond  the  aim  of  our  work  and  therefore 
it  will  not  be  considered. 

The  heat  flow  provided  by  the  thermal  source  and  flowing 
through  the  ceramic  plate  is  [14]: 

C Ih=K(Th-Thj ),  (1) 

whereas  the  heat  flowing  from  the  thermocouples  to  the  external 
heat  sink  is: 

Qc  =  K(TCJ-Tc).  (2) 

The  voltage  generated  by  the  TEG  is: 

VT  =  N(  gcp  —  ocn)(THj  —  Tq)  =  aATeff,  (3) 

being  ap,  a„,  and  a  =  N(otp  -  an)  the  Seebeck  coefficients  of  p-  and  n- 
type  semiconductors  and  of  the  whole  TEG,  respectively,  and 
ocpn  =  (ocp  -  a„)  is  the  Seebeck  coefficient  of  the  single  thermocouple. 
In  case  of  open  electrical  circuit,  the  rate  of  the  heat  flowing  through 
the  thermocouples  is  simply  Q.in  =  Km(THj  -  Tq).  Let  us  consider 
now  a  closed  circuit  where  the  variable  load  resistance  RL  is  con¬ 
nected  to  the  TEG  electrodes  (see  Fig.  2)  and  consequently  an  elec¬ 
trical  current  7  flows  in  the  circuit.  In  this  case,  it  can  be  shown  [6] 
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Fig.  2.  Scheme  of  the  circuit  considered  in  the  physical  model.  VT=oiATeff,  Rin: 
internal  resistance,  RL:  variable  load  resistance. 
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Fig.  4.  Schematic  section  of  the  device  (box  +  heat  source  +  TEM). 


that  the  rate  of  heat  absorbed  and  the  rate  of  heat  dissipated  are 
respectively: 


Qh  =  ocTHjI  +  KjnATeff  —  -RinI2,  (4) 

Q .c  =  ttTql  +  KinATeff  +  - RinI 2,  (5) 


where: 


1  = 


Vt 

Rin  +  Rl 


(6) 


The  term  proportional  to  I  is  due  to  the  Peltier  effect,  whereas 
the  term  proportional  to  I2  takes  into  account  the  Joule  effect 
(the  factor  1/2  arises  from  the  boundary  conditions  of  the  heat 
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Fig.  3.  Picture  of  the  device  (box  +  heat  source  +  TEM).  (a)  Operative  shape,  (b) 
Open  view. 


transport  equation  [6]).  According  to  this  circuit  and  by  writing 
the  load  resistance  in  term  of  internal  resistance  as  RL  =  mRin  (m 
is  the  proportional  coefficient  [15])  the  output  power  PL  supplied 
by  the  TEG  is: 


PL  =  VLI  =  N2oc2ATi 


m 


Rin^+mY 


(7) 


where  VL  is  the  voltage  across  RL  and  Eq.  (6)  has  been  considered. 
Eq.  (7)  clearly  shows  that  the  maximum  power  is  obtained  when 
the  load  resistance  matches  the  internal  resistance,  that  is,  m  =  1, 
Rl  =  Rin- 

The  relationship  between  AT  and  ATefr  is  determined  by  consid¬ 
ering  that,  at  the  equilibrium,  the  heat  rates  flowing  through  the 
ceramic  plates  are  coincident  with  the  heat  rates  absorbed  and  dis¬ 
sipated  by  the  thermocouples,  giving  [14]: 


ATeff 


K 


I(  +  2  Kin  + 


2N2agnrM 

^in+^L 


AT  =  PAT, 


(8) 


where  Tm  ~  {Th  +  Tc)I 2  is  the  mean  temperature.  The  term  in  the 
denominator  containing  the  Seebeck  coefficient  describes  the  tem¬ 
perature  change  due  to  the  Peltier  effect,  which  depends  on  the  cur¬ 
rent  and  consequently  on  the  load  resistance.  However,  as  we  are 
showing  in  the  sequel,  this  term  is  small  and  can  be  reasonably 
neglected. 

The  TEG  efficiency  ri  is  defined  as  the  ratio  between  the  output 
power  PL  and  the  absorbed  heat  rate  QH  (see  Eqs.  (7)  and  (4)). 
Therefore  we  obtain: 


rl(m)=7T  =  -. 


2mZ(THj  -  Tq) 


Qh  2(1  +  m)  +Z(THJ  +  2  mTHJ  +  Tq) 
where  the  figure-of-merit  Z  is  defined  by: 

N2oc2 

_  n  (xpn 
KjnRin  ' 

The  efficiency  at  the  matched  load  (m  =  1)  is  equal  to: 


(9) 


(10) 


Table  1 

Parameters  of  the  thermoelectric  module  obtained  from  its  data  sheet. 


N  VDC 

Rin 

Imax 

Pmax 

Pmax 

W 

L 

H 

(V) 

(Q) 

(A) 

(W) 

(C) 

(mm) 

(mm) 

(mm) 

127  9.40 

4.30 

1.09 

5.10 

250 

40 

40 

3.8 

343 
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Fig.  5.  Scheme  of  heat  source  (top),  TEM  (middle),  and  heat  sink  (bottom).  THE: 
temperature  of  heat  source  at  the  catalytic  burner  side;  TH:  temperature  of  heat 
source  at  the  TEM  side;  TCe :  temperature  of  heat  sink  at  the  open  side;  Tc: 
temperature  of  heat  sink  at  the  TEM  side;  THJ :  temperature  of  hot  junction;  TCj : 
temperature  of  cold  junction.  I<H,  I<T,  and  I<D  represent  the  thermal  conductances  of 
metallic  link,  whole  TEM,  and  heat  sink,  respectively. 

_  1  -  Tq 

w  match  2  THj  -  \  (T„j  -  Tg)  +  \ ’  1  j 


a  slit  along  the  whole  basis  perimeter;  such  opening  ensures  the 
necessary  air  flow  required  by  the  heat  source,  without  loss  of 
functional  heat.  Two  additional  “wings”  have  been  fixed  on  the 
surfaces  of  short  faces  in  order  to  protect  the  TEM  and  the  heat 
source  from  external  unrequired  air  flow.  Inside  the  box  we  added 
polystyrene  to  stress  the  adiabatic  properties. 

3.2.  Heat  source 

As  heat  source  we  have  considered  a  flameless  catalytic  burner, 
which  in  our  prototype  is  a  simple  commercial  Zippo™  hand  war¬ 
mer  [18].  Such  burner  uses  fuel  for  lighters  which  releases  heat 
through  a  reaction  with  a  platinum  catalyst.  At  room  temperature 
and  in  an  open  space,  the  hand  warmer  cap  surface  reaches  the 
maximum  temperature  of  75  °C  ca.  The  hand  warmer  is  placed  in¬ 
side  the  box  with  the  burner  exiting  from  the  hole  of  the  upper  cap. 
In  the  box,  the  bottom  slit  and  the  upper  hole  guarantee  the  correct 
oxygen  supply  to  the  catalytic  burner. 

3.3.  Thermoelectric  module 


whereas  the  maximum  efficiency,  obtained  when 
m  =  mMAx  =  VI  +ZTm,  is  given  by: 


*7  MAX  — 


ThJ  -Tq  a/1  +ZTm  -  1 
Thj  VI  3~ZTm  +  j ^ 


(12) 


showing  that  the  Carnot  efficiency  rjc  =  (THJ  -  TCj)IThj  is  dropped  by 
a  term  which  depends  on  the  figure-of-merit. 


3.  Design  and  realization  of  TEG  device 

The  realized  device  is  characterized  by  the  following  parts:  (i) 
box;  (ii)  heat  source;  (iii)  TEM;  and  (iv)  external  energy  conversion 
circuit.  A  picture  of  the  whole  device  (without  external  circuit)  is 
shown  in  Fig.  3,  whereas  a  schematic  section  is  shown  in  Fig.  4. 

3.1.  Box 

The  aim  of  the  box  is  to  support  all  the  parts  of  the  device  and  to 
maximize  the  heat  transport  from  the  heat  source  to  the  TEM.  In 
particular,  it  is  a  handmade  rectangular  box  made  with  plywood 
sheets  1  cm  thick.  The  cross  section  is  rectangular 
(100  mm  x  98  mm),  coincident  with  the  cross  section  of  the  heat 
sink,  whereas  two  opposite  sides  are  longer  than  the  other  ones 
(100  mm  x  133  mm  and  98  mm  x  99  mm)  to  support  the  heat 
sink.  The  upper  cap  is  fixed  on  the  top  of  shorter  faces  and  it  has 
a  square  hole  in  the  middle  with  the  same  lateral  dimensions  of 
the  TEM.  The  lower  cap  is  fixed  to  the  box  body  in  order  to  leave 


The  TEM  used  in  our  device  is  the  commercial  GM250-1 27-14- 
16  of  European  Thermodynamics  [19].  Data  sheet  gives  the  param¬ 
eters  shown  in  Table  1.  TEM  is  tested  at  the  hot  side  temperature  of 
250  °C  and  at  the  cold  side  temperature  of  50  °C.  N  is  the  number  of 
thermocouples,  VDC  the  maximum  voltage  at  open  circuit,  Rin  the 
internal  resistance,  IMAX  the  maximum  current,  PMAX  the  maximum 
power,  and  TMAX  the  maximum  temperature  of  the  hot  surface.  W 
and  L  are  the  width  and  the  length  of  the  ceramic  plates,  respec¬ 
tively.  Direct  measurements  give  h  =  2.2  mm,  h!  =  0.8  mm,  and 
1=  1.5  mm. 

To  transfer  the  heat  from  the  catalytic  burner  to  the  hot  side  of 
the  TEM  we  used  an  aluminum  heat  sink  between  the  heat  source 
and  the  TEM.  The  heat  sink  is  fin  type  with  the  same  cross  section 
of  the  TEM’s  one.  The  fins  were  cut  with  an  appropriate  shape  to 
place  the  burner  close  to  the  TEM.  The  heat  rejected  by  the  TEM 
flows  through  a  heat  sink  placed  on  the  top  of  the  box,  fixed  to 
the  higher  box’s  faces.  This  heat  sink  constitutes  the  cold  side  of 
the  TEM.  Polystyrene  has  been  placed  around  the  TEM  to  thermally 
insulate  it. 

Let  us  discuss  how  to  optimize  the  relative  dimensions  of  hot 
side  heat  sink,  TEM,  and  cold  side  heat  sink.  By  considering 
Fig.  5,  we  introduce  the  temperatures  THE  (temperature  of  heat 
source  at  the  catalytic  burner  side)  and  TCe  (temperature  of  heat 
sink  at  the  open  side)  and  the  thermal  conductivities  I<H,  I<T,  and 
Kd  of  metallic  link,  whole  TEM,  and  heat  sink,  respectively.  In  par¬ 
ticular,  1  IKt  =  1  IK  +  1 1  I<in  +  1  IK  and  consequently  I<T  =  KinI<l 
(. I<  +  2 Kin).  By  assuming  a  perfect  thermal  insulation  along  the  lat¬ 
eral  faces  and  by  neglecting  thermal  convection  and  irradiation, 


Fig.  6.  Block  diagram  of  the  integrated  solution  with  LT  components. 
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in  open-circuit  condition  the  heat  conduction  along  the  TEG  is  gov¬ 
erned  by  the  following  equations: 


Kh(The-Th)  =  Kt(Th-Tc ),  (13) 

Kt(Th-Tc)=Kd(Tc-Tce).  (14) 

By  combining  the  two  equations  and  by  rearranging  the  terms, 
after  having  introduced  the  notations  A TE  =  THE  -  TCe  and  l//(eff  =  1  / 
I<H  +  1  II<T  +  1  II<D  and  the  additional  equation  I<D{TC  -  TCe )  =  I<eff {THE  - 
-  TCe ),  we  obtain  the  following  relationship: 


A  T  = 


1  +  I<H  +  I(D 


A  Te. 


(15) 


The  TEG’s  design  is  optimized  when  the  difference  between  AT 
and  ATe  is  as  small  as  possible.  Let  us  write  the  thermal  conduc¬ 
tance  as  I<i  =  kiAilhi,  i  =  H,  T,  D,  where  kit  Ait  and  h,  are  the  thermal 
conductivity,  the  cross  sectional  area,  and  the  height  of  the  corre¬ 
sponding  element  in  the  TEG,  where  the  shape  of  a  simple  rectan¬ 
gular  parallelepiped  is  considered  (for  sake  of  simplicity  in  the 
notation,  we  have  temporarily  introduced  a  subscript  for  the  geo¬ 
metrical  parameters  concerning  to  the  TEM).  Besides,  we  reason¬ 
ably  assume  that  kH  ~  kD.  Through  these  considerations  we  write: 


AT  =  Th  -  Tr  = 


1  i  h  ( At  ho  I  hA 
kH  hT  '  Ah  hT J 


A  Te. 


(16) 


A  further  consideration  should  be  done:  since  all  the  thermo¬ 
couples  have  to  receive  the  maximum  available  amount  of  heat 
without  any  dispersion  in  the  environment,  we  must  have  AT  ~  AH. 
Therefore,  for  given  kT  and  kH  we  have  to  fulfill  the  conditions  AT/ 
Ad  <c  1  and  hDlhT  <c  1  in  order  to  minimize  the  difference  between 
AT  and  A TE. 


Fig.  9.  Scheme  of  the  circuit  used  to  characterize  the  TEM.  VT=  aATeif,  Rin\  internal 
resistance,  RL  =  RT  +  fy:  variable  load  resistance.  R/  is  made  by  two  resistances  in 
parallel. 


3.4.  External  energy  conversion  circuit 

In  order  to  properly  manage  the  energy  produced  by  the  gener¬ 
ator,  it  is  necessary  to  consider  that  our  system,  when  switched  on, 
starts  to  produce  energy  till  the  end  of  the  fuel  also  if  no  load  uses 
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Fig.  8.  Prototype  of  the  energy  conversion  circuit. 
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Fig.  11.  Dependence  of  AT  on  the  load  resistance. 


Fig.  12.  Dependence  of  VL  (blue  dots,  left  axis)  and  I  (red  dots,  right  axis)  on  the  load 
resistance.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 


that  energy.  The  best  solution  is  not  only  to  convert  this  energy  to 
be  used  with  the  most  part  of  microprocessors  or  sensors  available 
on  the  market,  but  also  to  store  that  part  of  the  energy  that  is  not 
instantly  used.  Due  to  this  limitation  it  is  important  to  design  a 
control  system  that  can  monitor,  manage,  and  store  the  energy  in 
order  to  increase  the  efficiency  of  the  whole  system.  For  this  reason 
the  energy  converter  can  be  divided  into  three  different  subsys¬ 
tems:  (i)  DC/DC  regulator;  (ii)  control  system;  and  (iii)  storage  sys¬ 
tem.  By  using  thermoelectric  generators  with  a  low  temperature 
difference  between  the  two  sides  it  is  necessary  to  manage  very 
low  voltage  and  current.  Also  for  standard  DC/DC  boost  converter 
this  is  not  easy  due  to  the  very  low  amount  of  voltage  and  energy 
available.  Thus,  we  have  to  increase  low  voltages  (<500  mV)  above 
the  minimum  threshold  required  by  the  boost  converter  (>2V). 
Since  our  thermoelectric  generator  produces  a  voltage  between 
50  mV  and  500  mV,  it  is  necessary  to  increase  more  than  40  times 
this  value  to  permit  a  good  conversion. 

In  order  to  simplify  the  construction  of  the  first  prototype,  an 
integrated  solution  is  preferable  to  reduce  cost  and  to  obtain  a  suit¬ 
able  industrial  solution.  For  these  reasons  some  solutions  from  Lin¬ 
ear  Technologies  (LT)  turned  out  to  be  useful  because  they 
integrate  into  a  single  chip  the  DC/DC  regulator,  the  charge  con¬ 
troller,  and  the  main  control  system,  by  reducing  cost  and  imple¬ 
mentation  time.  Only  few  additional  components  have  been 
selected  and  added  to  the  integrated  chip  in  order  to  obtain  the  fi¬ 
nal  required  solution.  An  integrated  circuit  that  permits  this  kind  of 
conversion  is  the  LTC3108  by  Linear  Technologies;  this  chip  is  a 
boost  converter  that  requires  only  some  external  components  to 
work  in  wide  conditions.  By  using  a  transformer  instead  of  a  gen¬ 
eric  inductance,  it  is  possible  to  work  with  very  low  input  voltages 
too.  In  particular,  our  prototype  is  equipped  with  coupled  induc¬ 
tors  with  a  turns  ratio  of  1:100  to  increase  the  input  voltage;  the 
second  stage  is  given  by  the  boost  converter  which  provides  an 
output  voltage  which  can  be  set  equal  to  2.35  V,  3.3  V,  4.1  V,  or 
5  V  to  properly  supply  any  kind  of  microprocessors  or  sensors. 
On  the  other  hand,  the  voltage  of  the  storage  system  in  fixed  at 
the  value  of  5.25  V.  Since  the  output  voltage  cannot  be  greater  than 
the  storage  voltage,  by  setting  the  output  voltage  at  the  lowest  va¬ 
lue  it  is  possible  to  maximize  the  available  stored  energy.  The  block 
diagram  of  the  whole  energy  conversion  circuit  and  the  scheme  of 
the  equivalent  electronic  circuits  are  shown  in  Figs.  6  and  7, 
respectively. 


Fig.  13.  Output  power  PL  normalized  by  AT2  as  a  function  of  load  resistance  RL.  Red 
points:  experimental  data;  blue  line:  data  fit.  Inset:  magnified  view  around 
maximum.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 


RL(Q) 


Fig.  14.  Efficiency  r\  as  a  function  of  load  resistance  RL.  Inset:  magnified  view 
around  maximum. 
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The  produced  energy  can  be  either  directly  used  by  an  external 
load  or  partially  stored  for  future  use.  The  storage  system  can  be 
divided  in  two  main  components,  based  on  different  technologies: 

•  Electrostatic  storage  (super  capacitors):  most  suitable  when  the 
external  load  requires  small  quantity  of  energy  in  short  inter¬ 
vals  (less  than  10  min)  and  the  generator  continuously  supplies 
the  average  energy  required.  The  auto-discharge  factor  is  higher 
than  in  the  chemical  solution  but  in  the  short  period  it  can  be 
neglected. 

•  Chemical  storage  (NiMH  batteries):  most  suitable  when  the 
generator  has  an  intermittent  production  or  the  external  load 
requires  medium/large  quantity  of  energy  with  long  pause 
interval  also  when  the  generator  is  switched  off.  The  auto-dis¬ 
charge  factor  is  negligible,  but  the  storage  efficiency  is  lower 
than  in  the  electrostatic  storage,  by  causing  efficiency  reduction 
during  the  charge-discharge  cycle. 

Both  solutions  have  been  tested  by  using  a  4.7  mF  16  V  and 
330  mF  5.5  V  double  layer  capacitors  for  the  electrostatic  storage 
and  a  battery  pack  of  150  mAh  4.8  V  for  the  chemical  storage; 
the  chemical  solution  becomes  useful  when  the  energy  absorbed 
by  the  load  is  characterized  by  a  long  interval  of  inactivity 
(>30  min).  However,  during  the  experiments  we  used  the  electro¬ 
static  solution  only,  by  connecting  the  4.7  mF  16  V  capacitor.  The 


picture  of  the  prototype  of  the  energy  conversion  circuit  is  shown 
in  Fig.  8. 

4.  Experimental  results 

Different  measurements  have  been  performed  to  physically 
characterize  the  TEM.  To  this  end,  one  has  to  apply  the  tempera¬ 
ture  difference  (TH  -  Tc)  and  subsequently  measure  the  output 
voltage  in  different  conditions.  To  correctly  realize  this  goal,  one 
should  set  TH  and  Tc  by  using  a  temperature  feedback  active  control 
in  order  to  keep  constant  both  these  two  temperatures  and  the 
mean  temperature  inside  the  TEM  during  the  specified  set  of  mea¬ 
surements,  to  reduce  the  effects  of  temperature  on  physical  quan¬ 
tities  [14].  However,  since  this  is  not  the  main  aspect  of  our  work, 
we  have  characterized  the  TEM  by  directly  using  the  prototype  of 
our  device,  where  both  TH  and  Tc  are  not  adjustable  to  a  fixed 
known  value.  On  the  other  hand,  we  prove  that  all  the  physical 
quantities  have  been  rigorously  obtained.  In  particular,  measure¬ 
ments  have  been  performed  by  considering  the  circuit  shown  in 
Fig.  9,  where  to  the  TEG  (composed  by  box,  heat  source,  and  TEM 
previously  discussed)  we  have  connected  a  passive  circuit  com¬ 
posed  by  a  variable  load  resistance  RT  and  a  fixed  resistance  Rh 
the  last  one  used  to  measure  the  current  through  the  correspond¬ 
ing  voltage  drop;  therefore,  the  previously  defined  resistance  RL 
corresponds  to  RL  =  RT  +  R h 


Fig.  15.  Dependence  on  time  of  measured  and  computed  physical  quantities  obtained  during  indoor  data  acquisition  at  room  temperature,  (a)  Tc,  TH,  and  TM’,  (b)  AT;  (c)  VL  and 
PL;  and  (d)  apn,  the  horizontal  blue  line  represents  the  value  of  Seebeck  coefficient  computed  from  power  fit.  (For  interpretation  of  the  references  to  color  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  16.  Dependence  on  time  of  measured  and  computed  physical  quantities  obtained  during  outdoor  data  acquisition  with  variable  temperature  and  environmental 
conditions,  (a)  Tc,  TH,  and  TM;  (b)  AT;  (c)  VL  and  PL;  and  (d)  ocpn,  the  horizontal  blue  line  represents  the  value  of  Seebeck  coefficient  computed  from  power  fit.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Computed  values  from  data  acquisition  are  so  AT,  TM,  I,  RL  =  VL/J, 
and  PL  =  VL I.  In  Figs.  10  and  1 1  we  show  the  variation  of  Tc,  TH,  and 
Tm  and  AT,  respectively,  recorded  during  the  experiment,  plotted 
with  respect  to  the  variation  of  the  load  resistance.  In  addition, 
in  Fig.  12  we  show  the  variation  of  VL  and  /  for  the  whole  RL  values 
considered  during  the  data  acquisition.  The  dependence  of  these 
temperatures  with  respect  to  RL  has  to  be  attributed  to  the  fact  that 
the  amount  of  heat  exchanged  depends  on  Peltier  and  Joule  effects 
too,  which  depend  on  the  load  resistance  since  they  are  functions 
of  the  current. 

By  plotting  PL  versus  RL  we  are  able  to  determine  the  required 
parameters  through  a  data  fit,  in  agreement  with  Eq.  (7).  On  the 
other  hand,  the  previously  discussed  equation  of  the  output  power 
as  a  function  of  the  load  resistance  has  been  obtained  by  assuming 
a  constant  temperature  difference  AT,  which  is  not  our  case  during 
the  whole  data  acquisition.  To  overcome  this  problem,  we  there¬ 
fore  plot  Pl/AT2  versus  RL  and  then  we  fit  the  values  with  the  func¬ 
tion  Pl/AT2  =  CiPL/(c2  +  Rl) 2,  being  the  fitting  parameters  equal  to 
Ci  =  ( Notpn/3) 2  and  c2  =  Pin,  where  we  have  taken  into  account  the 
relationship  ATeff=  PAT.  The  result  is  shown  in  Fig.  13. 

The  average  of  mean  temperature  inside  the  TEM  and  the  aver¬ 
age  of  temperature  difference  around  the  TEM  during  this  data 
acquisition  are  respectively  (TM)  =  320.06  K  and  (AT)  =  1 1.35  K 
(see  Figs.  10  and  11,  respectively).  Data  fit  provides 
Ci  =  1.6631  mW  Q/I<2  and  c2  =  Rin  =  4.18  Q,  which  is  in  good 
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Fig.  17.  Voltages  measured  during  the  test  of  the  whole  device  as  functions  of  time. 
Blue  line:  VL  multiplied  by  10;  green  line:  Vout;  red  line:  Vstore.  The  voltage  drop  and 
rise  at  3800  s  ca.  is  when  the  LED  has  been  connected  and  disconnected, 
respectively.  The  voltage  drop  at  6000  s  ca.  is  when  the  external  circuit  has  been 
disconnected  from  the  TEG.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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agreement  with  that  one  supplied  by  the  data  sheet.  In  particular, 
we  obtain  ps  =  16.82  pQ  m.  The  accuracy  of  this  fit  is  guaranteed 
by  the  chi-squared  test.  By  considering  an  error  of  0.001  W/I<2  for 
each  experimental  data  plotted  in  Fig.  13  (error  bars  are  not  visible 
because  inside  the  points),  we  obtain  as  reduced  chi-squared 
y2  =  0.47,  values  which  confirms,  being  y2  <  1,  the  accuracy  of 
the  fit  and  the  correctness  of  Eq.  (7). 

Let  us  now  discuss  how  we  can  determine  the  Seebeck  coeffi¬ 
cient  of  a  single  semiconductor  couple.  In  first  approximation,  p 
is  equal  to  (see  Eq.  (8)): 

'’“rrk-0-9834-  07) 

where  the  contribution  from  the  Peltier  effect  has  been  neglected. 
Therefore,  we  simply  have: 

Xpn=w= 326,52  ^V/Ki  (18) 

which  is  in  agreement  with  the  values  known  in  literature  (see,  for 
comparison,  [6,16,15]).  Conversely,  by  rigorously  applying  the  exact 
expression  of  p  given  by  Eq.  (8)  we  obtain  a  second  order  parametric 
equation  for  ocpn  with  RL  and  TM  as  parameters.  By  solving  such 
equation  for  each  data  set  we  obtain  several  values  of  apn,  varying 
from  327.64  pV/K  to  326.55  pV/K  and  with  a  mean  value  of 
326.96  pV/K.  Moreover,  the  corresponding  values  of  p  run  from 
0.9801  to  0.9833.  Since  the  difference  with  the  approximate  value 
is  irrelevant,  in  the  sequel  we  will  consider  the  simplified  expres¬ 
sion  for  p. 

Under  these  considerations,  the  computation  of  the  figure-of- 
merit  (Eq.  (10))  gives  rise  to  Z=  9.05  x  10-4 1/K,  whereas 


ZTm  =  0.29.  According  to  the  previous  discussion,  as  stated  by  Eq. 
(9),  the  efficiency  varies  along  the  range  of  load  resistance  values 
as  shown  in  Fig.  14.  In  details,  characteristic  values  are 
77match  =  0.2205%  and  ^max  =  0.2214%  (see  Eqs.  (11)  and  (12)) 
against  a  Carnot  efficiency  equal  to  y\c  =  3.43%,  since  the  averaged 
temperature  values  of  hot  and  cold  junctions  provide  Thj  ~Tm  +  ~ 
AT/2  =  325.64  K  and  TCJ  =  TM  -  AT/2  =  314.48  K.  It  is  clear  that  the 
efficiency  is  definitely  low,  but  this  fact  is  due  to  the  relatively 
low  values  of  temperature  difference.  With  the  same  figure-of- 
merit,  if  we  suppose  THJ  =  400  K  and  TCj  =  300  K  we  would  obtain 
^7  MAX  =  1.94%. 

After  this  first  set  of  measurements  developed  to  characterize 
the  TEG,  we  carried  out  a  different  type  of  data  acquisitions  under 
similar  operative  conditions  observed  during  the  applications  of 
the  thermoelectric  element  as  power  generator.  In  more  details, 
we  fixed  the  resistance  load  to  a  well  defined  value,  we  turn  on 
the  heat  source,  and  then  we  acquire  data  every  20  s  for  6  h.  Such 
kind  of  measurements  have  been  implemented  both  in  laboratory, 
at  room  temperature,  and  outdoor,  where  the  temperature  and  in 
general  the  environmental  conditions  are  arbitrary  and  obviously 
not  adjustable. 

The  results  obtained  in  the  first  case  are  shown  in  Fig.  15,  where 
all  the  physical  quantities  are  plotted  as  functions  of  time.  In  this 
case,  we  have  RL  =  6.66  Q.  After  a  fast  rise,  the  temperature  TH  at 
the  hot  side  progressively  but  slowly  decreases,  whereas  the  tem¬ 
perature  Tc  at  the  cold  side  remains  basically  constant  thank  to  a 
cooling  fan  close  to  the  heat  sink.  Consequently,  both  the  mean 
temperature  TM  and  the  temperature  difference  AT  have  the  same 
functional  behavior  of  TH  (see  Fig.  15a  and  b).  It  follows  that  also 
the  output  voltage  VL  and  the  output  power  PL,  after  a  fast  rise 


Fig.  18.  Operative  tests  in  environmental  conditions,  (a)  Value  of  Vout  (2.85  V)  in  open  circuit  condition;  (b)  LED  lightning;  and  (c)  environmental  site  of  operative  tests. 
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during  the  transient  initial  state,  decrease  as  expected  from  the 
behavior  of  AT  (see  Fig.  15c).  Besides,  the  corresponding  efficiency 
ri  from  a  value  close  to  0.22%  monotonically  drops  to  0.18%.  It  is 
important  to  notice  that  the  output  voltage  is  always  noticeably 
higher  than  200  mV,  a  critical  value  required  in  several  electronic 
applications,  as  discussed  in  the  sequel.  Finally,  from  the  relation¬ 
ship  VT  =  RinI  +  VL  =  NotPnPAT  we  were  able  to  compute  the  Seebeck 
coefficient  at  each  acquisition  time;  the  result  is  shown  in  Fig.  15d 
where  we  plotted  the  value  of  otpn  obtained  from  the  power  fit  pre¬ 
viously  discussed.  Basically,  the  variation  of  otpn  is  not  influent  and 
can  be  reasonably  attributed  to  a  change  in  the  measurement  con¬ 
ditions.  Such  dependence  on  time  of  all  these  quantities  can  be 
attributed  to  the  fact  that  the  hot  air  remains  close  to  the  output 
of  the  catalytic  burner  and  therefore  the  reduced  absorption  of 
oxygen  decreases  the  heat  production. 

The  same  type  of  measurements  has  been  carried  out  by  placing 
the  TEG  outdoor,  leaving  the  device  subjected  to  random  environ¬ 
mental  conditions  (alternation  sun-shadow,  wind,  humidity,  etc.). 
During  this  data  acquisition  we  set  RL  =  7.48  Q.  The  results  are 
shown  in  Fig.  16.  With  respect  to  the  previous  case,  here  we  cannot 
highlight  a  common  behavior  of  the  variation  of  the  physical  quan¬ 
tities  since  the  random  change  of  external  conditions  is  reflected  in 
a  random  variation  of  measured  and  computed  parameters.  In  par¬ 
ticular,  the  external  air  flow  supplies  oxygen  to  the  catalytic  burner 
by  increasing  the  efficiency  of  heat  production.  It  was  observed 
that  after  a  relatively  strong  air  flow  TH  immediately  decreases  be¬ 
cause  the  hot  air  is  removed,  but  later  TH  reaches  higher  values 
since  the  new  oxygen  increases  the  heat  production.  As  conse¬ 
quence,  all  the  other  parameters  reflect  the  same  behavior.  As 
one  can  appreciate  by  direct  comparison  of  Fig.  15  with  Fig.  16, 
the  performances  are  better  when  the  TEG  is  placed  outdoor,  in 
particular  the  efficiency  is  basically  higher.  Finally,  we  notice  that 
the  first  drop  at  1020  s  was  due  to  the  undesired  switch  off  of  the 
heat  source. 

5.  Operative  tests 

Laboratory  measurements  have  been  developed  to  test  the 
whole  device,  by  including  the  external  energy  conversion  circuit. 
The  results  are  shown  in  Fig.  1 7  where  we  measured  the  TEG  volt¬ 
age  VL  ( VTEG  in  Fig.  17  since  it  corresponds  to  the  generated  voltage 
in  open  circuit  conditions),  the  storage  device  voltage  Vstore  (in  this 
case,  the  4.7  mF  16  V  capacitor),  and  the  voltage  Vout  generated 
through  the  conversion  circuit,  set  at  2.3  V.  Once  the  catalytic  bur¬ 
ner  is  switched  on,  in  open  circuit  conditions  VL  goes  quickly  to  its 
maximum  value.  Then  we  connected  the  external  circuit:  immedi¬ 
ately  VL  decreases  and  Vout  reaches  the  set  value  since  it  supplies  a 
resistor  of  984  kQ,  whereas  Vstore  increases  up  to  5.25  V  ca.,  accord¬ 
ing  to  the  characteristic  time  of  the  capacitor.  Later,  at  3800  s  ca., 
we  connected  a  LED  (forward  voltage  2.1  V,  current  from  3  mA  to 
10  mA)  for  few  seconds  at  the  external  circuit:  in  that  period,  as  ex¬ 
pected,  both  Vout  and  Vstore  decrease  and  the  LED  is  switched  on  un¬ 
til  Vstore  >  Vout.  When  the  LED  is  removed,  both  Vout  and  Vstore  reach 
their  previous  values  according  to  their  characteristic  times.  Final¬ 
ly,  when  the  TEG  is  disconnected,  Vstore  starts  to  decrease  and  Vout 
quickly  drops  when  Vst0re  ^  Vout ,  whereas  VL  returns  at  its  original 
value.  Such  kind  of  experimental  test  proves  that  the  TEG  is  able  to 
cyclically  supply  low  power  devices,  represented  by  the  LED  during 
these  measurements,  with  a  duty  cycle  determined  by  the  charac¬ 
teristic  times  of  the  system. 

In  order  to  verify  the  capability  of  the  device  to  operate  in  real 
environmental  conditions,  we  tested  it  at  the  top  of  Testa  Grigia,  a 
summit  along  the  border  line  between  Italy  and  Switzerland  in  the 
Monte  Rosa  range  (coordinates:  45°56/2.03//N-7°42/29.94//E; 


altitude  above  the  sea  level:  3480  m;  date:  May  25,  2012  -  see 
Fig.  18).  Since  the  use  of  laboratory  instruments  was  not  easily 
available,  we  limited  the  experiment  by  measuring  Vout  with  a  por¬ 
table  multimeter  and  by  testing  the  LED  lightning.  Fig.  18a  and  b 
shows  the  value  of  Vout  (2.85  V)  in  open  circuit  condition  and  the 
LED  which  is  correctly  switched  on,  respectively.  Finally,  Fig.  18c 
shows  the  environmental  site  where  the  device  has  been  tested. 
To  our  knowledge,  this  is  the  first  case  in  which  this  type  of  device 
is  tested  in  real  operative  conditions,  ready  to  be  used  in  such 
hard-environmental  situations. 

6.  Conclusion 

In  this  work,  we  discussed  the  realization  and  the  characteriza¬ 
tion  of  a  thermoelectric  generator  prototype  used  to  power  auton¬ 
omous  sensors  placed  in  remote  environmental  locations.  A  first 
set  of  experimental  tests  have  been  carried  out  in  order  to  obtain 
the  physical  parameters  which  typically  characterize  the  thermo¬ 
electric  devices  (Seebeck  coefficient,  thermal  conductance,  internal 
electrical  resistance,  emitted  power,  etc.);  the  obtained  values  are 
in  agreement  with  those  ones  found  in  literature  for  such  devices. 
In  particular,  we  prove  that  the  physical  model  can  be  accurately 
checked  by  directly  making  use  of  our  prototype,  without  using  a 
further  experimental  setup  specifically  committed  to  this  aim.  Sub¬ 
sequently,  we  performed  a  series  of  measurements  directed  to  test 
the  capability  of  our  device  to  continuously  supply  energy  to  low- 
power  circuits;  the  results  show  that  the  prototype,  with  a  low- 
cost  catalytic  heat  source,  supplies  an  output  voltage  larger  than 
200  mV  and  an  average  power  close  to  9  mW  for  at  least  6  h,  when 
a  load  resistance  close  to  the  internal  resistance  is  connected  to  the 
thermoelectric  module.  Finally,  we  tested  the  whole  device,  given 
by  the  TEG  connected  with  an  appropriate  external  energy  conver¬ 
sion  circuit,  showing  that  it  is  able  to  supply  low  power  sensors.  In 
particular,  the  device  has  been  successfully  tested  in  real  hard 
environmental  conditions. 
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